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1. Introduction 

Supersymmetry is a key ingredient in the construction of models of fundamental physics, 
since it provides for a salient possibility to combine internal symmetries with the Poincare 
group. Even though distinguishing features of supersymmetric systems can be understood 
within perturbation theory, many important properties such as collective condensation 
phenomena often related to symmetry breaking are inherently nonperturbative. If super- 
symmetry is realized in nature, powerful and flexible nonperturbative tools will be needed 
to investigate the underlying mechanisms of these strong-coupling phenomena. 
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As supersymmetry does not only mix bosons and fermions but also involves spacetime 
translations, lattice methods built on spacetime discretization often go along with a partial 
loss of supersymmetry. The construction of appropriate lattice formulations in addition to 
the challenge of dealing with dynamical fermions is an ongoing effort ^, ^, Q|. These 
studies need to be complemented by nonperturbative continuum methods preferably with 
manifest supersymmetry. 

In recent years, the functional renormalization group (RG) has become such a non- 
perturbative tool as has been demonstrated by many successful applications ranging from 
critical phenomena, via fermionic systems and gauge theories even to quantum gravity, 
see Is], 0, |9|, for reviews. However, the number of applications to supersymmetric 
systems is rather small. In this work, we formulate and test the functional RG for a simple 
supersymmetric system, namely, supersymmetric quantum mechanics. 

In fact, ordinary quantum mechanics has often been used for illustrating and testing the 
nonperturbative capabilities of the functional RG, since the RG flow equations are easily 
obtained and approximate solutions can directly be compared to known exact results or 
high-precision numerics. In particular, the study of ground- and excited-state energies 
with RG techniques has received a great deal of interest in the last few years |11, 12, 13 



14| , y. Whereas single- well potentials can be treated comparatively easily even at extreme 



coupling, double-well potentials have turned out to be more challenging, since the analytic 
RG flow equations have to build up the non-analyticities from tunneling; the latter are 
usually described in terms of instantons, being of topological nature. 

In jl^, Horikoshi et al. study the quantum double well using an expansion in powers 
of the field, [^] and |14, [l^ go beyond this approximation and solve the RG flow in the so- 
called local-potential approximation for the effective potential (i.e., leading-order derivative 
expansion). Within the propertime RG, Zappala [^] also includes wave function renor- 
malization (i.e., next-to-leading-order derivative expansion), and finds good agreement for 
the mass gap. Particularly, this study convincingly demonstrates that the functional RG 
automatically includes also fluctuations of topological degrees of freedom without explicitly 
introducing them by hand. 

Supersymmetric quantum mechanics was introduced by Witten |16] as a toy model 
for spontaneous symmetry breaking. The first to use renormalization group methods for 
supersymmetric quantum mechanics were Horikoshi et al. |11|. They investigated a broken 
supersymmetric model with nonperturbative renormalization group methods and calcu- 
lated the nonvanishing ground-state energy and that of the first excited state in a poly- 
nomial expansion of the effective potential. They found good agreement with the exact 
results for all cases where tunneling is not important. This latter region has been cov- 
ered in |15] within the propertime RG, where again the observation was made that a wave 
function renormalization improves the results for the energy spectrum, i.e., helps including 
tunneling. 

Both approaches use regulators that break supersymmetry which makes it difficult to 
distinguish between explicit and spontaneous or dynamical supersymmetry breaking. One 
possibility to solve this problem is the inclusion of symmetry breaking by the regulator 
into the symmetry relations as done in [17, In fact even the lattice discretisation can 
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be viewed as a supersymmetry breaking regulator. The corresponding modified symmetry 
relations are similar to the Ginsparg- Wilson-relation, introduced in [19| and extended in 
|2C], and were established for super symmetric models in |21|. But so far a solution of these 
relations is possible only in some simple cases. In this paper, we present an approach to 
flow equations for supersymmetric quantum mechanics which maintains supersymmetry 
manifestly on the level of the RG flow equation with the aid of an invariant regulator. In 
contrast to [11| and [15|, we concentrate on a system with unbroken supersymmetry. 

Our approach is similar to the works by Bonini and Vian [^, 23] where a supersym- 



metric regulator for the id Wess-Zumino model is presented. The functional RG has also 
been formulated for supersymmetric Yang-Mills theory in |24| employing the superfield for- 
malism; for applications, see also |25, 26]. Very recently, Rosten has investigated general 
theories of a scalar superfield including the Wess-Zumino model with the aid of a Polchinski- 
type of RG equation with elegant applications in the context of non-renormalization the- 
orems ]|2^. A construction of a Wilsonian effective action for the Wess-Zumino model by 
perturbatively iterating the functional RG has been performed in ]28]. 

The paper is organized as follows: in Sect. ||, we briefly recall the basics of Euclidean 
supersymmetric quantum mechanics, also introducing our notation. In Sect. |3| we derive 
the flow equation for the superpotential and introduce a general class of supersymmetric 
regulator functions. In Sect. § we discuss the flow equation for the superpotential for 
different regulators. In Sect. ^ we introduce wave function renormalization and in Sect. ^ 
we compare our results with exactly known results. 



2. Euclidean supersymmetric quantum mechanics 

For our study of supersymmetric quantum mechanical RG flows, we employ the super- 
field formalism to maintain supersymmetry manifestly. The Euclidean superfield has the 
expansion 

^ = (j) + eiij + 4)9 + eeF (2.1) 

with anticommuting parameters 6, 9. Supersymmetric interaction terms are obtained as 
Z'-term of 

W{<^) = W((l)) + {9ip + ^9)W'{^) + 99{FW'i(l)) - W" {(p)i:il)) , (2.2) 

where the superpotential M^(<I>) is a polynomial in and W{(p) denotes the same polyno- 
mial evaluated at the scalar field cj). The nilpotent supercharges Q = idg + 9dr and Q = 
ide+9dr anticommute into the generator of (Euclidean) time-translations, QQ+QQ = 2idr- 
Supersymmetry variations are generated by = eQ — eQ, such that the variation of the 
superfield takes the form 

5,$ = e{i^p + i9F + 9^ + 9^) - [itp + i9F -9<j) + 99i)) e, (2.3) 

from which we read off the transformation rules for the component fields, 

6(j) = ieip - iipe , 6ip = {<j) - iF)e , 6ip = e{4> + iF), 6F = -eip - ipe . (2.4) 
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The super-covariant derivatives D = idg — Od-r and D = idg — 9dr fulfill similar anticom- 
mutation relations as the supercharges, 



{D,D} = {D,D} = {) and {D,D} 



-2idr 



(2.5) 



They commute with dr and anticommute with the supercharges. The integration over the 
anticommuting variables extracts the D-term of a superfield 



From this, we obtain the invariant action in the superfield formalism: 

^1 



(2.6) 



5[0,F,V',V'] 



drdOde 
1 



dr 



- i^p^ + iFW'{(l)) - i^pW"{(l))ip + -F^ 



(2.7) 



with kinetic operator K = \{DD — DD). Eliminating the auxiliary field -F, we obtain the 
on-shell action 



dr 



(2. 



It contains the bosonic potential V{(j)) = ^(H^'((/)))^ and a Yukawa term. In this paper, 
we consider models with unbroken supersymmetry. They have vanishing ground state 
energy = and are realized for superpotentials whose highest power is even. On the 
microscopic scale, we will focus on quartic superpotentials 

W{4^)=e4>+'^<P' + l<P^ + \<i>\ (2.9) 

as the defining starting point of the interactions of our quantum mechanical system before 
fluctuations are taken into account. 



3. Flow equation in the off-shell formulation 



3.1 Flow equation for the effective action 

The functional RG can be formulated in terms of a flow equation for the effective average 
action Tf, [^]. This is a scale-dependent action functional which interpolates between 
the microscopic or classical action S and the full quantum effective action T, being the 
generating functional for IPI Green's functions. The interpolation scale k denotes an 
infrared IR regulator scale which suppresses all fluctuations with momenta smaller than k. 
For A: ^ A with A denoting the microscopic scale, no fluctuations are included such that 
Ffc^A S. For A; ^ 0, all fluctuations are taken into account and we arrive at Tk~^Q — > F, 
i.e., the full solution of the quantum theory. The effective average action can be determined 



from the Wetterich equation |2£] 



dk^k 



r?^ + Rk 



dkRk 



(3.1) 
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which defines an RG flow trajectory in the space of action functionals with the classical 
action serving as initial condition. Here, F^^^ denotes the second functional derivative with 
respect to the dynamical fields, 

where the indices a, b in the general case summarize field components, internal and Lorentz 
indices, as well as spacetime or momentum coordinates. In the present case, we have 
= ((/>, F, -0, -f/^). (Note that ^' is not a superfield, but merely a collection of fields.) 
The supertrace in eq. (|3.lD as well as the pattern of functional differentiation in eq. ( ^I^ ) 
takes care of the minus signs from Grassmann- valued variables. The regulator function 
guarantees the IR suppression of modes below k, the shape of which is to some extent ar- 
bitrary; examples will be given below. Different 72^ correspond to different RG trajectories 
manifesting the RG scheme dependence, but the end point T^^o T remains invariant. 



The fiow equation (3.1) is an exact equation, involving the regularized exact propagator 



a one-loop structure. It can be viewed as the differential coun- 
terpart of a functional integral, or path integral in quantum mechanics. Its perturbative 
expansion yields full standard perturbation theory, but also nonperturbative systematic 
expansion schemes can be devised. In the present work, we use a derivative expansion of 
the effective action in powers of the covariant derivative in the off-shell formulation. This 
expansion is systematic in the sense that all possible operators can uniquely be classified, 
and it is consistent, since dropping higher-order terms leads to a closed set of equations. 
Most importantly, a truncation of such an expansion preserves supersymmetry. In this 
work, the derivative expansion of supersymmetric quantum mechanics will be worked out 
to next-to-leading-order. For simplicity, let us here begin with the leading order, corre- 
sponding to the local-potential approximation for the superpotential; to this order, the 
truncated effective action reads 



dr 



(3.3) 



The prime always denotes the derivative with respect to the bosonic field 4>. In the following 
we will derive flow equations for the superpotential Wfc(0). The next order which includes 
a wave function renormalization will be considered later on. 

Let us finally mention that the effective action is particularly convenient for extracting 
physical quantities: the effective action V = Tk=o evaluated on the solution of its quantum 
equation of motion yields the ground state energy, which is zero if supersymmetry is un- 
broken. Since T is the generating functional of IPI Green's functions, it provides access 
to all correlators and corresponding quantities. For instance, the location of the pole of 
the propagator T^'^\p'^) = is a measure for the energy of the first excited state in super- 
symmetric quantum mechanics (corresponding to particle masses in quantum field theory) . 
In the derivative expansion, this excited-state energy can directly be related to properties 
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of the superpotential, see below. An alternative to the effective-action flow would be the 
flow of the Wilson action Sk which has the advantage of being regulator-independent at 
leading-order in the derivative expansion [^], but can suffer from numerical instabilities 
within truncations [Q. 

3.2 Supersymmetric regulators 

For a supersymmetric initial condition and truncation, the flow and the resulting effective 
action is supersymmetric provided the regulator does not break the symmetry. When 



deriving the flow equation (3.1) from the functional integral, the regularization is introduced 



(2) 

by means of an additional action contribution ASk, such that = ■ The action 

principle therefore guarantees a supersymmetric regularization, as long as AS'fc is invariant. 
Indeed, an off-shell supersymmetric cutoff action can be written in terms of superfields and 
its covariant derivatives: 

^Sk = \j dT ^Rk{D, D)^-se . (3.4) 
Since D and D satisfy the anticommutation relations ( |2.5D the regulator can be written as 
Rk{D,D) = in{-dlk)+r2{-dl,k)K, K = ^{DD-DD). (3.5) 

The factor i in front of ri is chosen for convenience such that the corresponding cutoff 
action matches the mass term. Similarly r2 is chosen such that its cutoff action matches 
the kinetic term. Both functions are functions of —d^, i.e., of in momentum space. For 
this general class of regulators, the cutoff actions read 

ASk =\ j dTdOdO $ (in + r2K) = \ j £^'^(-p)i?fc(p)*(p) , (3.6) 

where = {<j), F,ijj,'ip). The quadratic form Rk{p) is block-diagonal, 




R, = I 1 with blocks Rl = ( ^"'^ '''\, Rl= ° + 



iri r2 I \P?"2 — tfi U 



(3.7) 



and hence does not mix bosonic and fermionic degrees of freedom. Three properties of 
the regulator Rk{p) are essential: (i) Rk{p)\p^ /k'^^o > in order to implement an IR 
regularization, (ii) Rk{p)\k2/p2~*o = 0) implying that the regulator vanishes for vanishing k, 
(iii) Rk{p)\k^A^oo oo which helps fixing the theory with the classical action in the UV. 

For manifestly supersymmetric cutoff actions ASk, supersymmetry relates the regula- 
tors of bosonic fields to that of the fermionic field. This puts further constraints on the 
admitted cutoff functions in a supersymmetric theory, as can be seen from the following 
example. In view of the regulator structure in eq. ( |3.7D , one may be tempted to set ri = 0. 
A natural choice for the regulator functions would then be such that the bosonic component 
~ p^r2 induces a gap for IR modes, e.g., r2(p^//c^) ~ fe^/p^ such that p'^r2 ~ k'^. Supersym- 
metry implies to the regulator pr2 for the fermions and to the regulator r2 for the auxiliary 
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field, both of which diverge in the IR for this choice. Even though regulators of this type 
are perfectly legitimate in the full flow equation, they lead to artificial IR divergencies at 
higher order in the derivative expansion, e.g., for a wave function renormalization. This 
problem can be avoided by a softer IR behavior of r2 and including a suitable nonvanishing 
ri. 

3.3 Regularized on-shell action 

The equation of motion for the auxiliary field in the presence of a cutoff is 

F = -^W\ = VK^(0)+ri0, h = l + r2, (3.8) 

where, for convenience, we have introduced the function h{p) and the shifted superpotential 
W containing the cutoff functions r2 and ri. The regularized non-local on-shell action 
becomes 



San = / dr 



-4>h<j> - iifjhtp - ii;W"((t)U + -W'((/>)iw'((^) 
2 2 n 

It is invariant under the following deformed supersymmetry transformations 



(3.9) 



5(t) = ie^-i^e, 5^ = [4)-'^W{(t))^e, # = e (^0 + i>V'(</>)^ . (3.10) 

These non-local transformations close on infinitesimal time translations, 

{5,,5,, - ,5,, ,5,,) (field) = 2i(e-ie2 - e-2ei)a,(field) , (3.11) 

provided the fermionic field satisfies the deformed Dirac equation hip + W = 0. With 
( ]3.9| ) we have constructed a regularized (nonlocal) on-shell action which is invariant under 
deformed supersymmetry transformations. 

Nevertheless, we would like to stress that the off-shell formulation is crucial for the 
construction of an invariant flow equation with one-loop structure. As the on-shell su- 
persymmetry transformations act nonlinear ly on the fields, the resulting cutoff action is 
not quadratic in the fields. Even though an on-shell supersymmetric flow can straight- 
forwardly be constructed from eq. (|3.9|), the resulting flow involves higher- loop terms and 
thus is much more difficult to deal with. 

3.4 Flow equation 

Returning to the off-shell formulation and using the block-diagonal structure of the regu- 



lator ( |3.7| ), the flow equation for the effective action Tk[4>, Fjipjip] written in component 
fields reads 

dkTk = ^ STr I [rf + Rk] OkRk^ = ^ Tr {dkRk Gk)BB " ^ Tr {d^Rk Gk)pp , (3.12) 

(2) 

where we have introduced the regularized full Green's function or propag ator Gk = (T):' + 



Rk) ^- Upon insertion of the truncation ( |3.3[ ) into eq. ( 3.12 ), we need to project only onto 
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the flow of the superpotential Wk- It can be done by extracting the flow of either the term 
Unear in F or the term proportional to tptp, cf. eq. ( |3.3D . This is a direct consequence of 
the manifest supersymmetry of this approach. As an illustration of this fact, we do it both 
ways. For the projection, it suffices to consider constant flelds, such that an expansion of 
the inverse Green's function in terms of the constant anticommuting spinors '0) "0 yields 

= +Rk^ Gq}^ + ^Mi + MsV' + V^Mg^/; . (3.13) 

The propagator itself reads 

Gk =Go,k - Go,k{i'Mi + M2tlj)Go,k 

+ Go,k {MiGo,kM2 - M2Go,kM, - M3) Go,ki^iP . (3.14) 

To proceed we use the block notation, 

N = f ""^A . (3.15) 

The nonvanishing blocks of the operators in the expansion (3.13) have the form 

/ V + iW'A I hp + iW"\ 

^W" h J' ^'^^^^^^^-[hp-tW" J' 

/yy\;(3) q\ I n iy\;(3)\ 

M^FB = -MiBF = Q J , M2BF = -M^fB = L q , (3-16) 



M^BB 







To calculate the full propagator Gk we must invert Gq \ . The inverse of Gq \ is block 
diagonal, and the diagonal blocks read for constant flelds 

{Go,u)bb = ^ (yj^,, ^^2 + ,^w(3)) and {Go,k)FF = ^ {G-,1)ff (3.17) 

with determinantal factors 

Af = /iV + {y^"? and Ab = Af + ihFW^^^ . (3.18) 



Since the regulator is block-diagonal, see (3^), only the diagonal blocks of the dressed 
propagator enter the flow equation ( p. 12 ). These blocks can be calculated with the help of 
( 3.14 ) . Inserting the regulator (^) flnally yields 

Str {dkRk Gk) = jdT {Ho + H3^PiJ) (3.19) 
with (p and F-dependent coefficient functions 

J 27r A_bA_p 
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and 



y 27r L^B^F 

The flow equation ( p.l2| ) relates the supertrace ( p.l9| ) to the variation of the effective action 
p.3| ). To project onto the flow for the superpotential, we differentiate the flow equation 
with respect to F and afterwards set F = = = 0. This yields 



dkWL 



2 dF 



r dp dkr2W - >v^'') + 2hdknW" 

F=o 2 J 2tt A| • I • J 



Integrating with respect to (and dropping the irrelevant constant of integration) finally 
yields the flow equation for the superpotential 

....... 1 fdp hdkn-dkr^W'icP) 

where we recall the abbreviations h = 1 + r2 and W" = ri + W^. This flow equation for 
the superpotential is one of the central results of our work. From the solution of ( 3.23| ), 



we can calculate the effective potential Vk by eliminating the auxiliary fleld in the effective 
action. In passing, we note that a quicker way to obtain the flow equation makes use of 
the superspace formulation, and an efficient approach is summarized in appendix |A[ 



The fiow equation ( 3.23 ) can alternatively be obtained by projecting the flow of the 



effective action ( 3.12 ) onto the coefficient of ijjtjj. This way one obtains 



dkW'^ = \H^\^^^. (3.24) 



The two projection formulas ( |3.22| ) and ( 3.24 ) indeed give rise to identical flows, since 



^\F=o = mF=o. (3.25) 
This identity illustrates the fact that our flow equation is manifestly supersymmetric. 

4. Flow of the superpotential for different regulators 

The regulator in the flow equation not only suppresses IR modes, but also guarantees 
UV regularization due to the operator insertion dkRk for Rk decreasing with momentum. 
This renders the flow local in momentum space, enhancing also the numerical stability. In 
quantum mechanics, this property is less important, since quantum mechanics is UV flnite. 
This allows to choose less UV-restrictive regulators for which the momentum integral in 
eq. (p. 23 ) can be carried out analytically. 



Indeed, as long as no diagrams with closed F loops contribute to the truncation, the 
regulator r2 can be dropped completely, as ri is sufficient to regularize all diagrams with 
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at least one (f> ov line, as is clear from the structure of the regulator (|3.7| ). Then the flow 
equation ( 3.23| ) simplifies to 



We verify in appendix that this regulator choice is sufficient for guaranteeing that 
the microscopic action is the correct starting point of the flow without closed F loops. 
Incidentally, setting ri = and using r2 as a regulator alone in the flow equations would 
lead to artificial divergencies for the wave function renormalization, as mentioned above. 
Next, we will discuss and compare different regulators. In principle, the choice of the 
regulator can be optimized in order to minimize truncation artifacts. However, due to 
the mixing between momentum- and field-dependencies in the denominator of eq. ( [4.1] ) 



~ ri{p^)Wl! {(/)), simple optimization strategies for bosonic systems [31, ^ do not apply 
and full functional optimization would be required However, since we are not aiming 
for high-precision calculations, our regulator choice will be guided by simplicity. 

4.1 The Callan-Symanzik regulator 

First, we consider a simple Callan-Symanzik regulator ri{p'^,k) = k for which eq. ( [4.1] ) 
reduces to the simple flow equation 

We will discuss and compare various approaches to solve this flow equation for different 
parameters and, in particular, for non-convex classical superpotentials. 

4.1.1 Polynomial expansion 

For a polynomial approximation, one may expand the superpotential Wk{(t>) in eq. ([4.2] ) in 
powers of the bosonic field <j), 

W,icP)=Y,^r with Wk^^i<P) = W,,i^) = e<P+'^<l>' + ^-<P^ + ^<p'. (4.3) 
^-^ n 2 3 4 

n 

Then also the right hand side of the flow equation can be expanded similarly. A comparison 
of coefficients leads to a system of coupled ordinary differential equations for the coefficients 
an{k). Terminating the expansions on both sides at order A'' and setting an>N the 
system becomes closed and can be solved numerically. At the cutoff k = A, the non- 
vanishing coefficients are (ai, a2, 03, 04) = {e,m,g,a). Note that for > 3ma the classical 
superpotential becomes non-convex. 

Indeed, such an expansion about = is not adjusted to the flow, as the largest 
contribution to the flow equation arises from field values which minimize W^'. An expansion 
of eq. ( [4.3[ ) about the minimum of W^!, 

W,{cP) = J2~^{^- Mk)T > W'^'i^o) = 2a3 = , (4.4) 

n=l 
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thus has a much better convergence behavior. At the cutoff, the initial conditions are 
provided by the nonvanishing parameters (oi, 02, 04, 0o) which can directly be linked with 
{e,m,g,a) given above. Most importantly, = 0,2 + 3d4{4> — 0o)^ is an even function 
of (j) — (j)Q. Thus, the flow is also even, implying that W'^ stays even at all scales and all 
coefficients an{k) vanish for odd n > 3. From the {(j) ~ 4>of coefficient of the flow, we find 
{k + a2)^(?A;0o = = which states that (jiQ is scale-invariant. The same is true for 

ai, since dtai = 0. The differential equations for the nontrivial even coefficients of the 
truncated system up to order = 10 read 



dka2 

dkae 
dkas 



3 a4 
9a| 



P 
SaeP 



k + a2 



p3 

3 27dl - SOdideP + ZagP^ 
2 P4 
81a| - ISBdldeP + {25dl + 420458)^^ 



p5 

5 243a| - 540a6a|P + [ISdasdj + 225dldi)P^ - TOagaeP^ - 5451004^"^ 
2 P6 



The energy Ei of the first excited state is determined by the curvature of the effective 
potential Vk = ^(VF^)^ at its minimum i^minj note that (f>min is generically not equal to 
(po- At the minimum, W' vanishes, such that Ei = V7"((/>min)- Table || contains the gap 
energy Ei for classical superpotentials with parameters e = m = a = I and different 
values of g. For > 3 the initial superpotential becomes non-convex. In addition. 



g 


0.0 


0.2 


0.4 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


1.8 




2.008 


1.960 


1.895 


1.815 


1.722 


1.615 


1.497 


1.371 


1.237 


1.097 




2.205 


2.140 


2.064 


1.980 


1.889 


1.794 


1.699 


1.608 


1.530 


1.472 




2.214 


2.146 


2.070 


1.987 


1.898 


1.808 


1.721 


1.646 


1.596 


1.590 




2.201 


2.135 


2.060 


1.977 


1.888 


1.798 


1.711 


1.638 


1.595 


1.612 


PDF 


2.203 


2.137 


2.062 


1.979 


1.890 


1.798 


1.710 


1.633 


1.584 


1.590 


exact 


2.022 


1.970 


1.905 


1.827 


1.738 


1.639 


1.534 


1.426 


1.323 


1.235 



Table 1: Energy Ei of the first excited state calculated in different orders of the polynomial 
approximation with the Callan-Symanzik regulator for e = m = a = 1. For comparison, also the 
results obtained from the solution to the partial differential equation ( |4.2| ) (PDE) and the exact 
values from numerically diagonalizing the Hamiltonian are given. 



the minimum (/>min moves away from the expansion point (po, in principle signaling the 
break down of the polynomial approximation which can be expected to hold only near 4>q. 
Nevertheless, the values Ei[g) obtained for the polynomial approximations of orders 4, 6, 8 
and 10 converge to values obtained by solving the full partial differential equation ([4.2|). 
We conclude that the polynomial expansion as an approximation to the full solution to 
leading-order derivative expansion works satisfactorily for the energy Ei at these coupling 
values. However, as the ~ 10% difference to the exact gap energies shows, the leading-order 
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derivative expansion itself gives acceptable but not very precise results. This should be 
compared to the analogous flow-equation approximation for non-supersymmetric quantum 
mechanics which yields an error below the percent level even at strong coupling. 

One important difference is that we have a flow equation for the superpotential and 
not for the effective potential itself. As a consequence, the flow equation tends to make 
the superpotential convex but not necessarily the effective potential. Figure |l| shows the 
flow of the effective potential Vk in the polynomial approximation (iA) with = 6 for a 
convex and non-convex Wc\- 



-0.6 -0.4 -0.2 




Figure 1: The effective potential 14 in the polynomial approximation for W^^{4)) — l + (j)+g(j)^ +4)'^ 
The left panel shows the potential for g = and the right panel for g = 2. 



4.1.2 Partial differential equation 

It is known from the study of non-supersymmetric systems that the polynomial approxima- 



tion fails for nonconvex potentials [12, 15]. The latter require a solution of the full partial 
differential equation ( [4. 2] ), which we did with NDSolve of Mathematica. In practice, we 
have chosen (f> in the range of G (—200, 200) and kept the potential at its classical values 
on the boundary of this range. The results for three different scales are depicted in Figure 
|2[ For convex superpotentials, the solutions obtained from the polynomial expansions and 
from solving the partial differential equation are almost identical. But in the non-convex 
case, the polynomial expansion fails to reproduce the correct asymptotic form of the super- 
potential. Non-convex classical superpotentials pose a numerical challenge as they might 
lead to instabilities originating from the singularity at W^{<j)) = —k. For such potentials - 
corresponding to a large coupling g - the flow equation also does not reproduce the correct 
gap energies Ei{g); see the PDF row in Table |l[ We shall see that similar conclusions hold 
for other regulators in the flow equation. 



4.2 Exponential and 6 regulator 

We want to compare the results obtained with the Callan-Symanzik regulator - which 
serves only as an IR regulator but does not suppress the UV ~ with an exponential and a 
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Figure 2: The effective potential obtained from tlic solution to the partial differential 
equation (4^) The left panel shows g ~ and the right panel g — 2. 



9 regulator, 



exponential regulator 



r^\p'^, k) = \/k'^ — • 9{k^ — P^) 9 regulator 



(4.5) 



In contrast to the infrared Callan-Symanzik regulator used in ([4.2|), these regularize the IR 
and UV. The corresponding flow equations for the superpotential read 



1 



dp 



c2 + 2/)e-pVfc^ 



-oo 

1 k 

4tt |/c2 - Ty^'2| 



vr (1 — sign Wk) + 2 arctan 



\k^-Wf\ 



(4.6) 



Note, that for the 9 regulator the integral (4.1) can be calculated analytically. The numer- 
ical results in Table |2| have been obtained from the solutions to these partial differential 
equations. For the exponential regulator we have taken cj) E (—20, 20) and the integration 
over p from —bk to bk. For the 9 regulator, we have used (p G (—50,50). The results 



g 


0.0 


0.2 


0.4 


0.6 


0.8 


1.0 


1.2 


1.4 


1.6 


1.8 


cs 


2.203 


2.137 


2.062 


1.979 


1.890 


1.798 


1.710 


1.633 


1.584 


1.590 


exp 


2.195 


2.130 


2.055 


1.972 


1.884 


1.791 


1.701 


1.622 


1.569 


1.684 


9 


2.197 


2.132 


2.058 


1.975 


1.888 


1.794 


1.705 


1.626 


1.576 


1.581 


exact 


2.022 


1.970 


1.905 


1.827 


1.738 


1.639 


1.534 


1.426 


1.323 


1.235 



Table 2: Energy of the first excited state for the classical superpotential (|4.3|) with (e, m, a) 



(1,1,1) and varying g calculated from the solution to the partial differential equation (4.1) with 
Callan-Symanzik, exponential and 9 regulators. 



for the three different regulators are depicted in Table S. They are almost identical, but 
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all differ on the ~ 10% level from the exact values displayed in the last row of the ta- 
ble. Higher precision thus requires a next-to-leading order calculation in the derivative 
expansion including a wave-function renormalization. 



5. Wave function renormalization 

To next-to-leading-order in the derivative expansion, a field-dependent wave function renor- 
malization is included in the truncation, 



,F,iP,'il;] = J drdede 
dr 



(5.1) 



with a field dependent function Zk{4>). The operator K has been defined in ( |3.5| ) and 
primes denote derivatives with respect to (j). The results of the last sections are recovered 
for = 

In the spirit of functional optimization we choose a spectrally adjusted regulator 
1 33, 34 1 which includes the wave function renormalization, 



ASk = ^J drdOde ^fc($)^^> (in + ra-fC) $ , (5.2) 



where Z'j^ is evaluated at a background field $ = {cp, 0, 0). The value of cj) can be viewed as 
a parameter labeling a class of regulator functions. In components, the cutoff action reads 

ASk = j dr Z'^,{4)f Q/c/tscA + ^FrsF + iFncj) + i){pr2 - iri)Tp^ . (5.3) 

Again, the flow of Zk can be read off from various operators. The simplest choice is given 
by the prefactor of the F"^ term, cf. eq. (^]^), since no time derivatives are involved here. 
After the projection onto the F"^ term at vanishing and a constant scalar field, we 
obtain the flow equations for the Callan-Symanzik regulator 



where we have introduced the abbreviations 

Af = {1 + kdk)Z'j^{^f and V = W" {(p) + kZl{$f . (5.5) 

To solve this system of coupled equations, we need to pick a value for the background 
field (j). Since we are interested in the excited-state energy, a reasonable choice would be 
4> = 4>mm- Since (pmin is not a priori known but a result of the flow, this would require an 
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iterative construction of the RG trajectory. Instead we make a technically much simpler 
choice and identify the background field (j) with the fluctuation field (p. Since all functions 
in the action are parameters of the background field (j), e.g., Zk{4>) ^ Zk{4>,4>), identifying 
(j) = 4> goes along with an approximation. This becomes obvious from the fact that, e.g., 
Z'kifpA = <P) = d<t,Zk{(t>,4>)\^=cf> / d^Zk{(t>,(t))- By setting = we ignore this latter 
difference. This approximation is well known in the context of background- field flows 
1 35, and the resulting flow can be viewed as a generalized propertime flow |37, |33| . 
As experience demonstrates, the error made by this approximation is outweighed by the 
improvement arising from the better spectral adjustment of the regulator, see, e.g., p8|| . 
Our results indeed confirm this conjecture. 

Including the wave function renormalization, the on-shell effective bosonic action at 
next-to-leading order in the derivative expansion is 



rfc[</.,V' = o,vi = o] 



dT 



Vk{4>) 



(5.6) 



At /c = 0, the energy gap results from the curvature of the effective potential with respect 
to canonically normalized fluctuations x = Z((p), for which we have the standard kinetic 
term ~ (c^rX)^- Hence, the energy of the first excited state for unbroken supersymmetry is 



lim 

fc->0 



ld^v,{Z'Hx)) 



:limi^ 

fc^O (^'(0))2 



(5.7) 



nin — ^{4^min) 

In Table |, the energy gap Ei{g) for (e, m, a) = (1, 1, 1) and various couplings g is compared 
with those obtained without wave function renormalization. The flow with wave function 
renormalization leads to much better results as compared to the flow without wave function 
renormalization. The agreement is very satisfactory with errors on the ~ 1% level even for 
couplings of order I. We conclude that the flow equation is able to capture nonperturbative 
physics in supersymmetric quantum systems with a reasonable precision. 
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1.584 


1.590 


PDE-hWF 


2.089 


2.031 


1.961 


1.879 


1.788 


1.690 


1.589 


1.489 


1.402 


1.341 


exact 


2.022 


1.970 


1.905 


1.827 


1.738 


1.639 


1.534 


1.426 


1.323 


1.235 



Table 3: Energy of the first excited state for the classical superpotential (4.2) with (e,m,a) — 
(1,1,1) and varying g calculated from the solution to flow equations with Callan-Symanzik regulator 
without and with wave function renormalization. 



6. Summary of the numerical results 
6.1 The energy of the first excited state 

We find that the polynomial approximation and the solution of the partial differential equa- 
tion without wave function renormalization for convex superpotentials converge to the same 
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Figure 3: We compare the energy gap Ei{g) computed with different approximation schemes for 
the classical superpotential Wd = f + mcf) + g^P' + (fi^ with to = f (left panel) and ra — i (right 
panel). For convex initial potentials, we obtain both a good convergence and a satisfactory accuracy 
of the next-to- leading order derivative expansion including a wave function renormalization. Beyond 
convex initial potentials, e.g., for larger couplings g > \/3 for m = 1, significant deviations from 
the exact result are observed, indicating a less controlled convergence behavior. 



value independent of the regulator, see Sect. [4.1.1| . Depending on the parameters of the 
classical superpotential, we obtain an accuracy of 10% for a small mass parameter (m ~ 1) 
and 2% for larger mass parameters (m ~ 3). Inclusion of the wave function renormahzation 
improves the results for the energy gap considerably. We achieve an accuracy of 3% for 
m = 1. Due to the presence of the auxiliary field, the wave function renormalization has 
contributions of order in the momentum and the F^-term - which is neglected without 
wave function renormalization - contributes to the on-shell potential 14(0). This effect is 
more pronounced for small mass parameters as the anomalous dimension scales with the 
inverse of m. For large m, the anomalous dimension is small so we do not expect large 
contributions in agreement with the numerical results. Figure |3| summarizes the results 
for the energies Ei{g) obtained from the different approximation schemes for m = 1 and 
m = 3. The explicit values are listed in Tables ||-^. 

The parameter space of large-a couplings is explored in Fig. ^ Here, we have used 
e = m = g = 1, implying that the initial potential is always convex. First, we observe that 
the excited-state energy from the polynomial expansion converges rapidly to that taken 
from the full solution at leading-order. The deviations from the exact result are again 
on the ~ 10% level. This is greatly improved at next-to-leading-order including the wave 
function renormalization. Here, the results match the exact values with an error on the 
1% level or below. The agreement holds over the whole coupling range from the weak- to 
the deeply nonperturbative strong-coupling regime. 

The overall picture confirms that the functional RG employing the super-covariant 
derivative expansion captures the physics of the first excited state well beyond the pertur- 
bative small-coupling regime. For initial boundary conditions given in terms of classical 
convex potentials, the derivative expansion appears to converge well and reaches a very 
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satisfactory accuracy level already at next-to-leading order. 

For combinations of couplings where 
the initial potential is non-convex, e.g., 
g > "v/S for e = m = a = 1, there is 
clearly room for improvements, as the 
deviations of the excited-state energy froi 
the exact result become large. Though 
the inclusion of a wave function renor- 
malization at next-to-leading order im- 
proves the result significantly, the accu- 
racy remams poor, see Fig. |. More- 
over, as the next-to-leading-order cor- 
rection becomes of the same order as 
the leading order, the convergence of the 
derivative expansion may become ques- 
tionable. On the other hand, it is im- 
portant to note in this context that the 
hierarchy of the derivative expansion is 
interwoven more strongly for the super- 
symmetric version than for non-super- 
symmetric systems. In the present case, 
also next-to-next-to-leading order operators can contribute to the flow of the superpoten- 
tial. These contributions may be relevant for non-convex initial potentials and thus restore 
the convergence properties of the derivative expansion. 

6.2 The global structure of the effective potential 

Whereas the polynomial expansion does rather well for the excited-state energy for the 
convex case, we observe its break-down beyond this restricted case: For instance for > 2> 
at e = m = a = 1, the classical superpotential ceases to be convex. Here, the polyno- 
mial approximation fails for asymptotic values of the field, since it tries to provide for a 
polynomial solution of the partial differential equation near the expansion point, where the 
low-energy effective potential Vfc becomes flat. The global structure of the effective poten- 
tial for g = 2 calculated from the partial differential equation (^]^) and the polynomial 
approximation with Callan-Symanzik regulator are plotted in Figure |5| together with the 
classical potential. 

As expected the polynomial approximation is not able to reproduce the correct global 
structure whereas the partial differential equation is able to do so. The other regulators 
lead to the same global structure of the effective potential. 

7. Conclusions 

In this paper, we have presented a functional RG approach to supersymmetric quantum 
mechanics. Our approach is formulated in terms of an exact and manifestly supersymmetric 
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Figure 4: Energy gap Ei{a) versus coupling a for 
e = m ^ g ~ 1 (convex initial potentials). We ob- 
serve a good convergence of the polynomial expan- 
sion. At next-to-leading-order derivative expansion 
including a wave function renormalization, satisfac- 
tory quantitative accuracy is obtained for the whole 
coupling range and even at strong coupling, demon- 
strating the nonperturbative capabilities of the func- 
tional RG. 




Figure 5: The effective potential W'{(f))'^ with the 
Callan-Symanzik regulator for a nonconvex Wd with 



flow equation for the effective action which is a supersymmetric variant of the Wetterich 
equation. We have used the supersymmetric off-shell formulation which is the crucial 
ingredient to maintain the simple one-loop structure of the flow equation. The approach can 
straightforwardly be generalized to other supersymmetric models based on a real superfield. 

We solve the flow equation nonper- 
turbatively in a systematic and consis- 
tent approximation scheme based on an 
expansion of the effective action in pow- 
ers of field operators of increasing num- 
bers of supercovariant derivatives. To 
leading order, this yields a flow equa- 
tion for the superpotential - a supersym- 
metric analogue of the local-potential ap- 
proximation; a fleld dependent wave func- 
tion renormalization appears in the flow 
to next-to-leading order. 

In the present work, we focus on un- 
broken supersymmetry by considering only 3 = 2- The polynomial expansion fails to repro- 

, , . 1 1 1 • u J. • duce the global structure of the effective potential, 

superpotentials whose highest power is 

whereas the full numerical solution of the superpo- 
even. As a physical observable, we con- ^ ^. , ^ „ ^ ^. 

tential flow agrees well with the expectation. 

centrate on the energy of the flrst excited 

state resulting from the effective potential. A comparison with the exact solution provides 
information about the convergence of the derivative expansion. Our results conflrm that 
the functional RG is indeed capable of describing the system over the whole range from 
weak to strong coupling. Our approach works particularly well for initial convex poten- 
tials. Here, first quantitative estimates can already be obtained from a simple polynomial 
expansion of the superpotential. For the excited-state energy, the polynomial expansion 
also converges nicely, whereas the solution of the full partial differential equation for the 
superpotential is required for global properties of the potential. Since the excited-state 
energy is a physical quantity, it should also be universal in an RG sense. In fact, our 
results show little dependence on the regulator which confirms this required universality. 
At next-to-leading order, the inclusion of a wave function renormalization improves the 
quantitative accuracy considerably. For convex potentials, the functional RG result agrees 
with the exact result within an error on the ~ 1% level even at strong coupling. 

As soon as the initial potential becomes non-convex, the fiow-equation result for the 
energy to lowest order starts to deviate significantly from the exact result. As is already 
known from standard quantum mechanics, the relevant tunneling processes are associated 
also with higher orders in the derivative expansion. Inclusion of the wave function renor- 
malization indeed improves our result, even though sizable deviations from the exact result 
still remain. The reason for this can be anticipated: supersymmetry forces us to organize 
the expansion in powers of the super-covariant derivative. This, however, mixes different 
orders of time derivatives; e.g, in the off-shell version of any supersymmetric theory with 
a scalar multiplet, the auxiliary field and the derivative of the scalar field occur on equal 
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footings. This is visible, for example, in the super symmetry transformation of ip being 
proportional to — iF, see (p.4|). On the other hand, we expect that the low-lying exci- 
tation energies are mainly determined by the long-wavelength fluctuations, such that an 
expansion in time derivatives of the field should be well justified. 

The crucial observation in this context is that the super-covariant derivative expansion 
contains terms without time derivatives also at higher super-covariant derivative order, for 
instance, ^[(DD)^]^ ~ -|- . . . . In particular, these F-potential terms can directly 
contribute to the flow of the superpotential. Since these terms are generated sizably only 
at larger values of the coupling, it is natural to expect that they can exert a pronounced 
influence on the energy gap at large coupling. As even higher-order operators will not take 
a direct influence on the flow of the superpotential, it is conceivable that the excited-state 
energy converges at this next-to-next-to-leading order of the super-covariant derivative 
expansion. Otherwise, the convergence and use of this expansion in the tunneling regime 
would be questionable. 

A study of these higher orders giving access to operators with higher powers of F 
are also needed for the case of broken supersymmetry. In this case, a nonzero vacuum 
expectation value of F is expected to occur, the description of which requires knowledge 
of the effective potential of this auxiliary field. 

The models considered here can be obtained by a dimensional reduction from the 2d 
Wess-Zumino model with A*" = 1 supersymmetry. This in part is the reason that most 
structural results of the present work also apply to this two-dimensional field theory, for 
example to the form of the cutoff action and the structure of the flow equations. The 
super-covariant derivative-expansion techniques are straightforwardly generalizable. Work 
in this direction is in progress. 
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A. The flow equations in superspace 

In this appendix we sketch the derivation of the flow equation for the superpotential in 
superspace. The equivalence of this manifestly supersymmetric derivation with the one in 
component form will be shown afterwards. The superspace-coordinates (x, 6, 9) are denoted 
by z. 

The supertrace that defines the flow of the effective action translates into a superspace 
integral: 

dkTk = \ I dzdz'duRk{z,z')Gk{z',z), = {vf + RkY^ (A.l) 
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As in the component formulation the fields are taken to be constant to calculate the Green's 
function Gk{z', z). In addition the expression is expanded in terms of the covariant deriva- 
tives D and D. To zeroth order in the covariant derivatives one finds 



J deS dtW{^) = \j^ deS de'S' {idtn{p) + dtr2{p)K{p)) x 



X sio' - me' - ^)^Tal^^(^"' - mo' - e) . (A.2) 

Note that in momentum space the operator K = \{DD — DD) still contains derivatives 
with respect to the Grassmann-coordinates. These derivatives act on the first entry of the 
adjacent delta- functions. The only two contributions that remain after an integration over 
9' and 9' are the ones where the highest Grassmann derivative acts on one and only one of 
the delta functions inside the integral. Therefore we get 



/ 



For the lowest component of the superfield this is exactly the fiow equation ( 3.23| ). 

To prove the equivalence of this derivation to the one in given in the main body of 
the paper we observe that the transition from component to superfield formulation can be 
achieved with the linear operators P{9, 9) = (1, 99, -9, 9) and P^((9, 9) = (1, 99, 9, -9) : 

^ = P,{9, 9){(l), F, {[,, V)^ = (1, 99, -9, 9U<P, F, i^, = cj) + 9^^ + i^9 + 99F 

= (</>, F, i,, i^Ul, 99, 9, -9)i = (</., F, ^, ^)iPT{9, 9) . (A.4) 

In the other direction the operator Q{9, 9) = {99, 1, —9, —9) must be applied: 

(0, F, i), = j d9d9 Qi{9, m = j d9d9 {99, 1, -9, -9)i {(p + 9^; + i'9 + 99F) . (A.5) 

Note that as expected P0' ,9')Qi{9 ,9) = Qi{9',9')P^{9,9) = 5{9' - 9)5{9' - 9) and 
/ d9d9Qi{9,9)Pj{9,9) = f d9d9 Pf {9,9)Qj{9,9) = 6ij. The operator Rk and its inverse 
can be easily translated from component to superspace formulation using these operators: 

{Rk{x,x'))ij = j d9d9d9'd9' P'[{9,9)Rk{z,z')Pj{9',9') (A.6) 

{Rk{x,x'))r^ = [ d9d9d9'd9' Q,{9,9){Rk{z,z'))-^Qj{9' ,9') , (A.7) 



with / dz' Rk{z, z'){Rk{z' , z)) = 5{z — z'). So the flow equations translate into 
dxdx'{-lY'dtRk{x, x')ijGk{x', x)ji 

d9d9 d9'd9' dxdz"{-lf^Qi{9, 9)Pi{9' , 9'){dtRk){x, 9, 9; z")Gk{z"; x, 9', 9') 
dzdz"{dtRk){z,z")Gk{z";z) (A.^ 
with (—1)'^' —1 if i is a fermionic index and 1 otherwise, since {—1Y'-Pi{9,9) = P^{9,9). 
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B. Initial conditions 



Throughout this work, we have set the regulator component r2 = 0. With regard to the 
regulator structure ( p.T] ) , one may wonder whether this choice is compatible with a sufficient 
UV suppression of all modes. If not, the initial condition of the flow would not necessarily 
coincide with the microscopic (classical) action, but a separate UV renormalization would 
be necessary. 

Indeed, it is easy to see that diagrams containing closed F loops with a momentum- 
independent free propagator can give rise to UV divergencies signaling this insufficient UV 
suppression. On the other hand, closed F loops do simply not contribute to the present 
truncation; this would require, e.g., the occurrence of F self-interactions ~ F^ which are 
generated only at higher-order in the super-covariant derivative expansion. Perturbatively, 
they occur at the two-loop level. We conclude that there is no danger from F loops up to 
next-to-leading order in the derivative expansion. 

Indeed, sufficient UV suppression can directly be verified. For large A;, the cutoff action 



AS'fc dominates the action in the defining Euclidean path integral which is of the form |29| 



(B.l) 

The integral becomes dominated by small fluctuations around the classical solutions in the 
presence of the cutoff. A good estimate is thus provided by a saddle-point approximation 
of the path integral. Using the simple Callan-Symanzik regulator ri = A; as an example, 
one-loop corrections are given by 

00 „ oo 

1 /■ dp , S^4,Sff -S],. I f dp^ / iFW" 

1 I -n ' / In I I I 

"Y _ C* _ 



Rescaling p with k yields 



oo 

fe,iloop = -k / — In 1 + 



27r V k'^p'^ + {W"/k + iy 



-oo 



This integral vanishes for A; — > oo so that no UV counterterms are necessary to define the 
initial conditions. The starting point of the fiow equation is indeed the classical action. 
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